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Abstract 

We calculate NLO QCD corrections to the lightest neutral Higgs boson production associated 
with top quark pair at hadron colliders in the minimal supersymmetric standard model(MSSM). 
Our calculation shows that the total QCD correction significantly reduces its dependence on the 
renormalization/factorization scale. The relative correction from the SUSY QCD part approaches 
to be a constant, if either Ms or nig is heavy enough. The corrections are generally moderate(in 
the range of few percent to 20%) and under control in most of the SUSY parameter space. The 
relative correction is obviously related to rrig, At and /i, but not very sensitive to tan/3, Ms at both 
the Tevatron and the LHC with our specified parameters. 



PACS: 12.60.Jv, 14.80. Cp, 14.65.Ha 



I. Introduction 



One of the major objectives of future high-energy experiments is to search for scalar Higgs particles and 
investigate the symmetry breaking mechanism of the electroweak interactions. In the standard model 
(SM) 1 , one doublet of complex scalar fields is introduced to spontaneously break the symmetry, 
leading to a single neutral Higgs boson h°. But there exists the problem of the quadratically divergent 
contributions to the corrections to the Higgs boson mass. This is the so-called naturalness problem. 
One of the hopeful methods, which can solve this problem, is the supersymmetric (SUSY) extension 
to the SM. In these extension models, the quadratic divergences of the Higgs boson mass can be 
cancelled by loop diagrams involving the supersymmetric partners of the SM particles exactly. The 
most attractive and simplest supersymmetric extension of the SM is the minimal supersymmetric 
standard model (MSSM) [21 E]. In this model, there are two Higgs doublets H\ and H2 to give masses 
to up- and down-type fermions. The Higgs sector consists of three neutral Higgs bosons, one CP-odd 
particle (^4°), two CP-even particles (h° and H°), and a pair of charged Higgs bosons (H^). 

However, these Higgs bosons haven't been directly explored experimentally until now. The pub- 
lished experimental lower mass bounds for the Higgs bosons presented by LEP experiments are: 
M h o > 114.4 GeV (at 95% CL) for the SM Higgs boson, and for the MSSM bosons M h o > 91.0 GeV 
and M40 > 91.9 GeV (at 95% CL, 0.5 < tan/3 < 2.4 excluded). The SM fits to precision electroweak 
data 4 indirectly set a limitation of the light Higgs boson, M^o < 200 GeV, while there should has a 
scalar Higgs boson lighter than about 130 GeV in MSSM. This lightest Higgs boson with mediate 
mass is certainly in the exploring mass range of the present and future colliders, such as the Tevatron 
Run II, LHC and LC. At a LC the cross section for e + e~ — > tih is small, about 1 fb for s/s = 500 GeV 
and rrih = 100 GeV [HI El- But it has a distinctive experimental signature and can potentially be 
used to measure the top quark Yukawa coupling in the intermediate Higgs mass region at a LC with 
very high luminosity. S. Dawson and L. Reina calculated the NLO QCD corrections to e + e~ —* tih 
process at LC's in Ref. |S]. And in references [51 1101 ITT] the SM electroweak corrections to the process 
e + e~ — > tih are calculated. H. Chen et al., have studied the QCD and electroweak corrections to 
the process 77 — > tih in the SM at LC's|12|. All these works show that the evaluation of radiative 
corrections is a crucial task for all accurate measurements of tih production process. 

There are various channels which can be exploited to search for the Higgs boson h° with intermedi- 
ate mass at TeV energy scale hadron colliders, such as gluon-gluon fusion Higgs boson production(<?g — > 
h°), the associated production with a weak intermediate boson (qq' — > Wh° , Zh°). Recently, the 
production channels pp/pp — > tih + X attracted the physicist's attentions, because these chan- 
nels offer a spectacular signature (W^+H^&B&B)^! an d provides a possibility in probing the Yukawa 
coupling |14H 15j. The total cross section for pp/pp — > tih + X at tree level and NLO QCD corrections 
in the SM have been studied in Refs.[H El III EEZ1 EH- 

The supersymmetric (SUSY) electroweak corrections to the e + e~ — » tih process can be over ten 
percent for favorable parameter values [TU]. In Ref.|2U], it was found that the SUSY QCD interactions 
by exchanging gluinos and squarks can impact on the Yukawa coupling vertex in the process e + e _ — ► 
tih at LC. At pp/pp hadron colliders with a center-of-mass energy of TeV scale, the dominated 
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contributions to tih° production are from subprocesses qq,gg — > tih°. To these high energy tth 
production processes, the SUSY radiative corrections, especially the SUSY QCD corrections, may be 
remarkable. 

In this paper, we calculated the cross section for the associated production of the Higgs boson 
with top quark pair in the MSSM at hadron colliders including the NLO QCD corrections. In section 
2, we present the calculations of the leading order cross sections to pp/pp — > tth + X in the MSSM. 
In section 3, we present the calculations of the 0(af) QCD corrections to pp/pp — ► tth + X in the 
MSSM. The numerical results and discussions are presented in section 4. Finally, a short summary is 
given. 



II. The leading order cross sections 

The Feynman diagrams at leading order(LO) for the subprocess 

q(pi)q(p2) - t(h)t(k 2 )h°(k 3 ), (2.1) 

in the MSSM are plotted in Fig.l. They are s-channel, gluon exchange diagrams with Higgs boson 
radiation off top-quark and anti-top-quark, respectively. The process 

9(pi)g(P2) - i(fci)t(fc 2 )/i°(A;3), ( 2 - 2 ) 

at the tree level in the MSSM are described by the Feynman diagrams of Fig. 2. The LO Feynman 
diagrams for both subprocesses in the MSSM are the same with their corresponding ones in the SM. 
In above two channels we use p\^ and fci 2,3 to represent the four-momenta of the incoming partons 
and the outgoing particles, respectively. Because of the small mass of u- and d-quark, we neglect the 
diagrams which involve h° — u — u and h° — d — d Yukawa vertexes. 




h° 

h°\ ' 

Fig. 1. The tree- level Feynman diagrams for the qq — > tth subprocess. 
The explicit expression for the amplitudes of subprocess qq — > tth at tree level can be written as: 

Mf Q = ^ fc (p 2 ) 7 ^ i ( Pl )^ % (fc 1 ) ^^j^ 7^Kfc2)^^ + (h - k 2 ) (2.3) 

where s = (pi +P2) 2 , T a is the SU(3) color matrix, A = ggY^ USY \ g s is the strong coupling constant. 
Y t [^ USY ^ is the Yukawa coupling between Higgs boson and top quarks in the MSSM. As we know the 
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Fig. 2. The tree-level Feynman diagrams for the gg — > tth subprocess. 
h° — t — i Yukawa coupling in the SM Y^ M ^ is expressed by 



(SM) 



tth 



2my/ 



But in the MSSM, y£? USY) is given as 



(SUSY) 



tth 



-i>9v 



mt cosa 
2mw sin/3 ' 



(2.4) 



(2.5) 



where a is the mixing angle which leads to the physical Higgs boson eigenstates h° and H°. The angle 
(3 is defined as tan/3 = V2/V1, where v\ and t>2 are the vacuum expectation values. 

According to the different topologies of Feynman diagrams, the explicit expression for the ampli- 
tudes of subprocess gg — > tth in the MSSM at tree level can be divided into three parts. 



M 99 _ is" 1 + M" 2 4- /If" 3 
■'"tree _ m tree + ■'"tree ^ lvl tree 



(2.6) 



where Mf^ e , Mf® 2 e and Mf^ e correspond to the amplitudes of Fig.2(a-b), Fig.2(c-f) and Fig.2(g-h), 
respectively. For the amplitude parts M^ e (i = 1,2,3), we have the expressions as: 



Mf r t = Af^T^ih^^g^ + (pa - pi) V " 2p&" A l 



fe 2 -h + m t 



Vi(k 2 ) + (h <-► k 2 ) (2.7) 



(2.8) 
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M tree = -l AT ik T kj U j( k vfa 7 7 fl v i( k 2) + (Pi P2, Cl £2) (2.9) 

— fci • P2 — «2 - Pi 

where and €2 are the polarization four- vectors of the incoming gluons. The SU(3) structure constants 
are given by / a & c . Then the lowest order cross sections for the subprocesses qq, gg — > tt/r in the MSSM 
are obtained by using the following formula: 



a 



11,99 
LO 



L- 1 d$ 3 ElCf I 2 (2.10) 



2|fci 

where d<&3 is the three-body phase space element. The summation is taken over the spins and colors 
of initial and final states, and the bar over the summation recalls averaging over the spins and colors 
of initial partons. The LO total cross section of pp/pp — ► tih° + X can be expressed as: 

a LO (pp/pp^tth° + X) = Y,J dx 1 dx 2 G p t (x 1 ,fi)G p j /p (x 2 ,fi)a i l (x 1 ,x 2 ,fi) , (2.11) 

ij 

where Sr % i {ij = qq,gg) is the LO parton- level total cross section for incoming i and j partons, Gf^'s 
are the LO parton distribution functions (PDF) with parton i in a proton/antiproton. 

From the above deduction, we can see that the ratio between the tree level cross sections of 
subprocess qq(gg) — > tih° in the SUSY model and the SM, is written as 



^Lo SY \qq,99 -> tthP) = cos 2 a 
*{ S o M \qq,gg^tih°) "sin 2 /?' 



(2.12) 



III. NLO QCD Corrections in the MSSM 

In the calculation of the NLO QCD corrections in the MSSM, we adopt the dimensional regularization 
in D = 4 — 2e dimensions to isolate the ultraviolet (UV), infrared(IR) and collinear singularities. 
Renormalization and factorization are performed in the modified minimal substraction(MS') scheme, 
and the wave functions of the external fields, and top quark's mass in propagators and in the Yukawa 
couplings are renormalized in the on-shell(OS) scheme. We divide the 0{o? s ) QCD correction to the 
subprocess qq(gg) — > tih° in the MSSM into two parts. One is the so-called SM-like QCD correction 
part, another is SUSY-QCD correction part arising from virtual gluino/squark exchange contributions. 
Then the total NLO QCD corrections and relative corrections in the MSSM can be expressed as 

A ~ (11,99) _ a - (11,99) , a - dim) 1(11,99) - xdi,99) , $(11,99) (o 1 \ 

NLO ~ SM-like ' SQC'D' ~ SM-like + °SQCD- 

where we define the relative correction as 6 = . The NLO SM-like QCD correction part (relative 

correction part) in the MSSM has following relation with the NLO SM QCD one 

COS CX 2 A A 

SM-like = (-^~o) Aa S M, ^SM-like = &SM- (3.2) 
sm jj 
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In our calculation we introduce the following counterterms. 



m t 



t L 



UL 
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mt + Smt, 
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1 + -5Z U L ) u L , 



9s + 5g s 
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1 + 



di, d R 



1 + ) u R 

1 + ^ ] 



(1 + -<5Z S )G M , 



(3.3) 



where <7 S denotes the strong coupling constant, t, u, d and denote the fields of top-, up-, down- 
quark and gluon. The definitions and the explicit expressions of these renormalization constants can 
be found in Ref. [U- For the renomalization of the QCD coupling constant g s , we use the MS scheme 
except that the divergences associated with the colored SUSY particle loops are subtracted at zero 
momentum (22]- Since we have 5g = dg( SM - llke ) _|_ §g( s Q CD ) ; the terms should be obtained as 



Sg 



(SM-like) 



An 



4' 



{SM-like) 



1 1 , mi 
- + - In — | 

e 3 jii*. 



(3.4) 
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where 
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An 
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^ 12 
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In 



(SQCD) 



>-i(n/ + l), 



,(3.5) 



(3.6) 



N = 3 ,nf = 5 and 1/e = 1/ejjv ~ IE + hi(47r). The summation is taken over the indexes of squark 
and generation. The MS 1 scheme violates supersymmetry explicitly, and the qqg Yukawa coupling g s , 
which should be the same with the qqg gauge coupling g s in the supersymmetry, takes a finite shift 
at one-loop order as shown in Eq. ()3.7|) 23 . 

n 4 

g s =g s [l + -l(-N-C F )}, (3.7) 

3 and Cf = 4/3. In our numerical calculation we take this shift between g s and g s into 



with N -- 
account. 

Actually, the calculation of the NLO SM-like QCD corrections in the MSSM for the subprocesses 
qq,gg — ► tih is the same as that of the NLO SM QCD corrections in Refs.[TH I15j . except their 
numerical results satisfied the relations shown in Eq.(3.2). 
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The NLO SUSY-QCD contribution part to the qq(gg) — ► tih° subprocess comes from the one- 
loop diagrams involving virtual gluino/squark exchange. For demonstration, we show the pentagon 
diagrams which contribute to the NLO SUSY-QCD correction part for the subprocesses qq — ► tth° 
and gg — ► in Fig.3, where the upper indexes s,t,u = 1, 2. Because there is no massless particle in 
the loop, all these diagrams with gluino/squark loop are IR finite. The pentagon and box diagrams 
in SUSY-QCD part are UV finite, but the self-energy and vertex diagrams in this part contain UV 
divergences. That is renormalized by the proper related counterterms defined in Eq.(3.3). 




Fig. 3. The pentagon diagrams for the qq —>■ tth° and gg — > tth° subprocess. 

The 0{o%) supersymmetric QCD correction part of the cross section in the MSSM to the subpro- 
cesses qq, gg — > tih° can be expressed as 

A4S = ^| / dd> 3 £2Re (M^M^) , (3.8) 
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where M^f are the Born amplitudes for qq,gg — ► tth° subprocesses, and M$qcx) are the renor- 
malized amplitudes of all the one- loop Feynman diagrams involving virtual gluino/squark. 

In the calculations of loop diagrams we adopt the definitions of one-loop integral functions of 
Ref . |24| . The Feynman diagrams and the relevant amplitudes are generated by FeynArts 3(22], and 
the Feynman amplitudes are subsequently reduced by FormCalc32. The phase space integration is 
implemented by using Monte Carlo technique. The numerical calculations of integral functions are 
implemented by using developed LoopTools. 

We write the NLO QCD corrected parton-level total cross section a^ LO {xi, X2, At) as: 

±.{m,gg) _ ~(qq,gg) , \±.{gg,gg) 
u nlo — °lo ^"nlo ' 

The NLO QCD corrected total cross section of pp/pp — ► tih° + X in the MSSM can be expressed 



as: 



ctnlo(pp/pp -> tth°) = J dx 1 dx2G p i (x 1 ,fi)G p j /p (x2,fJ,)a ( ^ i l (x 1 ,x 2 ,fJ,) , (3.9) 

where &x{} (ij = qq,gg) is the NLO QCD corrected parton-level total cross section for incoming i 
and j partons, and G V J P are the NLO parton distribution functions (PDF) for parton i in a pro- 
ton/antiproton. The equation include two channels: qq,gg — ► tihP. In our calculations, we choose the 
factorization scale equals the renormalization scale, i.e., fj,t = fx r = Q. The partonic center-of-mass 
energy squared, s, is given in terms of the total hadronic center-of-mass energy squared s = X1X2S. 



IV. Numerical Results and Discussion 

In our numerical calculation, we adopt the MRST NLO parton distribution function |28j and the 2-loop 
evolution of a s (/j?) to evaluate the hadronic NLO QCD corrected cross sections, while for the hadronic 
LO cross sections we use the MRST LO parton distribution function and the one-loop evolution of 
a s (fi 2 ). We take the SM parameters as cw^l)" 1 = 127.918, m w = 80.423 GeV, m z = 91.188 GeV, 
wit = 174.3 GeV , m u = rrid = 66 MeF[2nj. There we use the effective values of the light quark masses 
{m u and m^) which can reproduce the hadron contribution to the shift in the fine structure constant 
a ew (m^) 30 . The other relevant parameters, such as mixing angle of the Higgs fields a and masses 
of the lightest Higgs boson, gluino, stop-quarks, are obtained by adopting the FormCalc package, 
except otherwise stated. The input parameters for the FromCalc program are Ms, M2, A t , m^o, n 
and tan [5. The related parameters for the MSSM Higgs sector are obtained from the CP-odd mass 
m^o and tan/3 with the constraint tan/3 > 2.5. In the program the grand unification theory(GUT) 
relation M\ = (5/3) tan 2 9\yM2 is adapted for simplification and the gluino mass rrig is evaluated by 
rrig = a s (Q) /a ew (mz) sin 2 9w M2. For the sfermion sector, the relevant input parameters are Ms, Af 
and fjL, and there we take the assumptions of Mq = Mjj = Mb = Me = Mi = Ms and the soft 
trilinear couplings for sfermions q and I being equal, i.e., A q = A\ = Af. 
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We present the dependence of the cross section on the renormalization/factorization scale Q/Qq 
in Fig.4(a-b) for the Tevatron and the LHC separately, where we denote Qq = rat + and the 

input parameters are taken as A t = 800 GeV, M s = 400 GeV, M 2 = 110 GeV, m A o = 270 GeV, 
H = —200 GeV and tan/3 = 6. With these input parameters, we get all the other supersymmetric 
parameters, among them cos a = 0.954, m h o = 120 GeV, ra^ = 207.75 GeV and = 577.63 GeV, 
but the value of m g is a function of the energy scale Q(ra g {Qo) = 317.9 GeV). In order to show 
the cross section dependence on the renormalization/factorization scale, we fix m g = 300 GeV in 
Fig.4(a-b). There we plot the curves for cross sections <tlo, &nlo and cr^J l4fa of the processes 
pp/pp —* tih° + X. The notations onlo an d c^xo'** 6 represent the cross sections involving complete 
QCD and SM-like QCD corrections. Fig. 4(a) shows that the NLO QCD contributions to the process 
pp —> tih° + X in the MSSM at the Tevatron, in which the dominant subprocess is qq — ► tih°, has 
a negative NLO QCD corrections near the position of Q = Qq. While Fig. 4(b) shows that the NLO 
QCD contributions to the process pp — ► tihP + X in the MSSM at the LHC, in which the dominant 
subprocess is gg — > tth°, will give positive corrections near the position of Q = Qq. Here we should 
note that if Q goes down to a very low value, i.e., Q « Qq, large logarithmic corrections spoil the 
convergence of perturbation theory in the proton-antiproton colliding energy of the Tevatron. That 
can be seen from our numerical results for the Tevatron. It shows that the total NLO QCD corrected 
cross section ctnlo i n the MSSM tends to have a negative value when Q — > 0|T3J . From Fig.4(a-b), 
we can conclude that the dependence of the NLO QCD corrected cross section ctnlo ° n the scale Q 
is significantly reduced comparing with ctlo, and is slightly weakened comparing with crff^Q hke . 

In the following calculation, we fixed the value of the renormalization/factorization scale being Qq. 
In Fig. 5(a) and (b) we show the LO and total NLO QCD cross sections ulo and ctnlo i n the MSSM 
as the functions of tan/3(m/jo) at the Tevatron and the LHC respectively, taking At = 800 GeV, 
M s = 400 GeV, M 2 = 110 GeV, m A o = 270 GeV and = -200 GeV. The corresponding relative 
corrections 5 of both cross sections versus tan/3(m^o), where the relative correction is defined as 
S = aNL °~J LO i are plotted in Fig. 5(c). From these figures, we can see that the cross sections ctnlo 
and ctlo decrease rapidly as tan (3 varies in the range from 2 to 10, and then goes down very slowly 
when tan/3 changes from 10 to 40. We can read from Fig.5(a-b) that when tan/3 increases from 2 to 
40, the total NLO QCD corrected cross section gnlo i n the MSSM decreases roughly from 9.1 fb and 
1078 fb to 5.2 fb and 641 fb for the Tevatron and the LHC, respectively. The two curves of relative 
corrections 5 for the Tevatron and the LHC in Fig. 5(c) look like rather stable when tan/3 runs from 
2 to 50. We can read from this figure that the NLO QCD relative correction values in the MSSM 
at the Tevatron and the LHC are generally about —17% and 26% in these varying range of tan/3, 
respectively. 

In Fig. 6, we show the relative NLO QCD correction S in the MSSM as a function of Ms, taking 
A t = 800 GeV, M 2 = 110 GeV, m A o = 270 GeV, fi = -200 GeV and tan^ = 6. The figure 
demonstrates that the relative NLO QCD corrections in the MSSM at the Tevatron and the LHC, are 
stable when Ms changes from 400 GeV to 2 TeV. Their values are about 25% at the LHC and —18% 
at the Tevatron. We find from our calculation that when Ms is taken as a large value, the correction 
from the NLO SUSY QCD correction part decreases to a constant due to the decouple of heavy stop 
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quarks. 

Fig.7 shows the total QCD relative correction 5 in the MSSM as a function of m g with the input 
parameters same as in Fig.4. We can see from this figure that the S have a concave structure in the 
vicinity of m g ~ 140 — 150 GeV , where the masses satisfy the relation nig + m t ~ + m h o and 
the re-scattering enhancement t\ — > 5 + i — > ii + /i° takes place. When m g goes from 400 GeV to 
2000 GeV, the relative corrections are very stable, they are about 24% for the LHC and —18% for the 
Tevatron. Similar with the case in Fig.6, due to the decouple effect the correction of the SUSY QCD 
correction part decreases to a constant when g is getting heavy. 

Fig. 8 presents the total NLO QCD relative correction 5 in the MSSM as a function of the SUSY 
parameter A t , assuming M s = 400 GeV, M 2 = 110 GeV, m A o = 270 GeV, fi = -200 GeV and 
tan/? = 6. We can see from the figure that the total NLO QCD relative corrections in the MSSM 
are very sensitive to A t in the region near the position of A t = 1000 GeU(where = 108.6 GeV). 
Actually, the reason for that is because of the large contribution from the light stop quark t\ loops. 
When the chosen parameters At and n, make a large mass splitting between the two scalar top-quarks, 
then the t\ becomes light. We can read from the figure the total NLO QCD relative correction 5 in 
the MSSM can reach -24% at the Tevatron and 7% at the LHC when A t is near 1000 GeV. 

In Fig.9, we show the total NLO QCD relative correction 5 in the MSSM as a function of the 
SUSY parameter p, assuming A t = 800 GeV, M s = 400 GeV, M 2 = 110 GeV, m A o = 270 GeV and 
tan (3 = 6. We can see that the total NLO QCD relative corrections in the MSSM increase slowly 
when jj, goes up from —1000 GeV to 1000 GeV , this is because the absolute values of the negative 
corrections from the SUSY QCD part are becoming smaller as /x increases. The value of 5 at the 
Tevatron can be beyond —22% when \i is about —1000 GeV . 

In this paper we calculated the NLO QCD corrections to the processes pp/pp — > tth° + X in the 
MSSM at the Tevatron and the LHC. We analyzed the dependence of the corrected cross sections or 
relative corrections on the renormalization/factorization scale Q, SUSY parameters tan/3, Ms, m g , At 
and [/,, respectively. It shows that the dependence of the total NLO QCD corrected cross section in the 
MSSM on the scale Q is significantly reduced comparing with the ulo- With our chosen parameters, 
the numerical results demonstrate that the relative correction is obviously related to m g , A t and \i in 
some parameter regions, but not very sensitive to tan (5, Mg at both the Tevatron and the LHC for 
our specified parameters. We conclude that the total NLO QCD corrections are generally moderate, 
which have the values in the range of few percent to about 20% in most of the SUSY parameter space. 
We find that the relative correction from the NLO SUSY QCD correction part becomes to be constant 
when either Ms or m g has large value. We find also the relative correction of the SUSY QCD part 
will be largely enhanced when the mass splitting between stop-quarks is large, the total NLO QCD 
relative correction in the MSSM 5 can reach —24% at the Tevatron and 7% at the LHC. 

Acknowledgments This work was supported in part by the National Natural Science Foundation 
of China and special fund sponsored by China Academy of Science. 



10 




Fig. 4. The cross sections aio a t the leading order and <jnlo involving the NLO QCD corrections 
in the MSSM as the functions of the renormalization/factorization scale Q with rrig = 300 GeV and 
the other parameters are from FormCalc by using the input SUSY parameters: Ms = 400 GeV, 
M 2 = 110 GeV, A t = 800 GeV, m A o = 270 GeV, n = -200 GeV and tan/? = 6. Fig.4(a) is for the 
process pp — ► tth° + X at the Tevatron and Fig.4(b) for the process pp — ► tt/t° + X at the LHC. 
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Fig. 5. The cross sections olo at the leading order and cjnlo involving the NLO QCD corrections 
in the MSSM as the functions of tan/3 with the input parameters At = 800 GeV, M2 = 110 GeV, 
m A o = 270 GeV, M s = 400 GeV, n = -200 GeV. Fig.5(a) is for the process pp -► tih° + X at the 
Tevatron and Fig. 5(b) for the process pp — > tth° + X at the LHC. Fig.5(c) shows the relative NLO 
QCD correction as the function of tan/3 in both the Tevatron and LHC. 
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